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Abstract 

We investigate the vacuum pp-wave and Aichelburg-Sexl-type solutions in f{R) and the modified 
Gauss-Bonnet theories of gravity with both minimal and nonminimal couplings between matter and 
geometry. In each case, we obtain the necessary condition for the theory to admit the solution and 
examine it for several specific models. We show that the wave profiles are the same or proportional 
to the general relativistic one. 
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I. INTRODUCTION 



Extensions of the general theory of relativity in which functions of some geometric quan- 
tity are coupled either minimally or nonminimally to the matter part of the usual action for 
gravity have been widely studied recently, mainly as part of current efforts to explain the 
Universe with its observed late-time accelerated expansion [l-3|. The so-called f{R) gravity 
comprises the largest subset of models constructed in this way, see Refs. j4-6| for reviews. In 
these models, a priori arbitrary functions of the scalar curvature of spacetime are included 
in the gravitational action in its various metric. Palatini, or metric-affine formulations. The 
f{R) function can also be coupled nonminimally to the matter Lagrangian, as suggested 
in Refs [3, B]. An interesting consequence of such nonminimal coupling is the emergence 
of extra forces making the trajectories of otherwise-free particles nongeodesic jo]. Another 
well-known class consists of the so-called J^{Q) gravity models, in which a function of the 
Gauss-Bonnet invariant, Q, is added to the Einstein-Hilbert action [l^. More general mod- 
ified Gauss-Bonnet theories of gravity with nonminimal coupling have also been suggested 
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[12], see also Ref. [l3|. 

Various aspects of the above-mentioned models have been studied in recent years, includ- 
ing black hole solutions and their thermodynamics and cosmological solutions with accel- 
erated expansion. There is now a huge literature on this, and a selection is listed in Refs. 
J:6| and IJ]- Massless and massive gravitational wave solutions have been presented, too 



151]. The issue of linearized f{R) gravity has been studied in Refs. 16Nl8|. Compared with 
general relativity, several new features arise in this context, namely, the appearance of extra 



19| 



polarization modes and a nonlinear dispersion law, see, e.g., Ref. 

Different models introduced in this context have been examined against several theoretical 
or observational criteria. As an example, the study by Dolgov and Kawasaki 2o| has shown 
that certain f{R) models accommodate ghosts and hence are ruled out by the consequent 
instability. One can also mention constraints coming from solar system effects such as 



planetary orbits or bending of light 
conditions 22|, in this regard. 



2l[ |. and the bounds coming from imposing the energy 



The aim of the present work is to find classes of these extended theories of gravity 
which admit plane-fronted-parallel rays gravitational wave solutions the same or similar to 
the general relativistic vacuum or Aichelburg-Sexl solutions 23|]. One motivation behind 
this work is the recent interest in gravitational pp-waves in general, and in gravitational 



shock waves of various sources, in particular. In fact, in the context of general relativity, the 



problem has been studied by boosting the Kerr metric in Refs. 
^err-Newman metric in Ref. 



3, 



2J] and [25|, by boosting the 



26|, in the presence of a nonvanishing cosmological constant 
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for 



for motion in Schwarzschild-Nordstrom and Schwarzschild-de Sitter spacetimes 
particles with arbitrary multipoles in Ref. 29|, for motion in Nariai universe 30|], and for 
motion in the presence of electromagnetic fields in Ref. [sjj. A reason for interest in such 



solutions lies in their role in the scattering of particles off of each other at ultrahigh energies 
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33| . Extensions of the Aichelburg-Sexl solution outside the context of general relativity 



3^ and 



have also been of interest, namely, in the framework of the brane- induced gravity 
in the context of the ghost-free model of massive gravity |35|. 

Second, and at the same time, this study aims to put forward a new criterion for testing 
extended gravity models. The basic idea is that the gravitational radiation by moving 
particles, massless ones in the present case, is detectable, at least in principle, and whether 
a given model admits the plane wave solution distinguishes it from other models. For models 
admitting the solution, the explicit form of the wave profile provides further information. 

In the next sections, after a brief review of the Aichelburg-Sexl solution of general rela- 
tivity, we apply the above ideas to several well-known extended theories of gravity including 
minimal and nonminimal f{R) gravities and J-'{Q) gravity. 



II. THE AICHELBURG-SEXL SOLUTION 



The Aichelburg-Sexl (AS) solution, first introduced in Ref. [23|, represents the gravita- 
tional field of a massless particle moving in an otherwise empty space. It belongs to the 
general class of plane-fronted gravitational waves with parallel rays, or pp-wave, solutions. 
The general form of a pp-wave line element may be written as 

ds = —dudv — K{u, X, y)du + dx + dy , (1) 

in which u = t — z and v = t + z are the null coordinates, and K{u, x, y) is obtained from 
the field equations. For vacuum, the field equation results in the two-dimensional Laplace 
equation to be satisfied by K{u, x, y). 

The AS solution may be obtained by inserting the energy- momentum tensor of a massless 
particle in the Einstein equation. Starting with the action for the coupled gravity-massless 



particle system, 



one obtains 



where 



^=7^ I V^Rd'x+ I ^C,(fx, (2) 



G^u = SttGT^i,, (3) 



is the scalar curvature, G^i, is the Einstein tensor, and 

rj. 2_6{Cpy^ 

is the energy-momentum tensor. The above action is expressed in a system of units in which 
c = 1. For convenience, we also set SttG = 1. For a massless particle of momentum p with 

'^P = ^ j 9^yX^'x''6^{x - x{T))dT (4) 

in which x^ = we have 

T"" =p j 5\x - x{T))x^'x''dT, (5) 

where x^{t) corresponds to the trajectory of the particle which has to be a null geodesic of the 
spacetime under consideration. Thus, we can take x^(r) to be of the form (r, 0, 0, 2;(r) = r) 
which is null and satisfies the geodesic equation in the spacetime described by Eq. ([1]). 
Inserting these together with the ansatz given by Eq. ([T]) into Eq. ([3]) and making use of 
the following relation 
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+ —]\n^/x^ + y^ = 2n5{x)5{y), 



we obtain 



P 



KAs{u,x,y) = 5(M)lnA/x2 + 2/2 (6) 



TT 



which upon substitution in Eq. ([T]), describes a gravitational shock- wave propagating with 
the speed of light along the 2;-direction. 



III. GRAVITATIONAL WAVES IN NONMINIMAL f{R) GRAVITY 
We start with the following action {q] 

s = j ^Tg + (1 + \F{R))C^ d'x (7) 

in which f{R) and F{R) are arbitrary functions of R. This represents an extension of 
the general theory of relativity in which the matter field Lagrangian density Cm is coupled 
nonminimally with the geometric structure F{R), and the coupling constant A controls the 
strength of the coupling. This reduces to the usual f{R) theories of gravity for A = 0, which 
in turn reduces to general relativity by choosing f{R) = R. 
The field equation associated with the above action reads 

S^u = V,V,f\R) - g,.Df'{R) 

+2A(V^V. - g^,n)CmF'{R) 

-2\F\R)CmR^, + (1 + XF{R))T^,, (8) 

where V means covariant differentiation, □ = V^V^, R^u is the Ricci tensor, and S^i, = 
f'{R)R,, - \f{R)g,.. 

Now, let us examine the above field equation to see if the vacuum pp-wave solution is 
admitted in this model. For vacuum, the last three terms in the right-hand side of Eq. ([8]) 
vanish, and if we insert the ansatz ([T]) into the resulting equation, we will reach the following 
relations 

/'(O) ( 7^ + 7^ ) ^(^' ^' y) - fiO)K{u, X, y) = 0, (9) 



and 



/(O) = 0. (10) 



Thus, the vacuum pp-wave solution is admitted only if the above condition holds. This 
condition is in fact among the requirements needed to constrain the models from cosmo- 
logical or solar system tests, which is also consistent with the limiting case of the ACDM 
phenomenology [36]. The other requirement is that f{R) should tend to a constant when 
the scalar curvature tends to infinity. 

When Eq. (fTOl) holds, Eq. ([9]) reduces to the two-dimensional Laplace equation in the 
transverse plane, provided /'(O) ^ 0. This means that the waveform is the same as the 
vacuum general relativistic one. For the special case where /'(O) = 0, Eq. ([9]) is trivially 



satisfied with arbitrary K{u, x, y). But this corresponds to the absence of the hnear term R 
in the Lagrangian and hence is ruled out. 

Now, we consider a massless particle moving along the 2;-direction and seek solutions of 
the form given in Eq. ([T]) (renaming K{u, x, y) to KnR{u, x, y) to avoid confusion). Inserting 
Eqs. (jl]) and (jSD into the above field equation (with replaced by Cp) and noting that for 
massless particles g^ux'^x'^ vanishes, we reach 



VlK^Riu,x,y) = -2pi±^^(5(«)5(x)5(y), (11) 

provided Eq. holds and /'(O) 7^ 0. Here, V| = ^ + Thus, if Eq. ^ is satisfied, 
the pp-wave solution ([1]) is admitted, and the waveform is given by 

KnR{u, X, y) = ^ ^^^^^^^ Kas{u, X, y). (12) 

If we also choose i^(0) = 0, then the effect of the nonminimal coupling disappears totally. 
For, F{fi) 7^ the nonminimal coupling has a contribution equal to jr^^- Because of the 
very small expected value of the coupling constant A, this would be a small contribution; 
see Ref. jsT] for a discussion of bounds on the values of A. 

It should be noted here that for the pp-wave spacetime ([T]), the nonlinear field equation 
(IH]) reduces to a linear equation discussed above, and this allows the use of distributional 
expressions as in general relativity. 



IV. GRAVITATIONAL WAVES IN f{R) GRAVITY 

The well-studied f{R) theories of gravity are in fact a subclass of the nonminimal theory 
considered above, with A = 0. Thus, both the vacuum pp-wave and the AS solutions are 
admitted if Eq. ^ holds. 

To determine the waveform of the AS solution, we insert A = into Eq. (fT2l) (this time 
with Kr{u, X, y) in place of KnR{u, x, y)). This yields 

KR{u,x,y) = j^KAs{u,x,y). (13) 

Since in general f'{R) > 0, otherwise ghosts are allowed; the signs of K{u,x,y) and 
Kas{u, X, y) are the same. 

An example of the models satisfying the requirement given by Eq. ( fTOl) is the broken 
power-law model 

c i—Y 

f{R) = R-m^ 



suggested in Ref. [36\, in which ci,C2,m, and n are constants, and n > 0. Note that the 
hnear term above is included to reproduce the Einstein- Hilbert action in Eq. ([7]). For this 
model, we have 



1 



lim , , 
R^o f'{R) 



1 if n > 1 
if n = 1. 



l-ci 

if n < 1 



By inserting this into Eq. f|T3|) . we conclude that for n > 1, the model admits the plane 
wave solution, with a wave profile the same as the general relativistic one for n > 1, and 
times the general relativistic waveform for n = 1. For n < 1, in which the above limit 



equals zero, the solution reduces to the Minkowski spacetime. In other words, for n < 1, 
the plane wave solution is not admitted, and this is in agreement with what we expect from 
the general requirement f'{R) > mentioned earlier. 

Another example is the Starobinsky model 38| given by 



Here, we have 



1 



/'(O) 



1, 



and hence the solution is exactly the same as the general relativistic one. 
Also, the following cosmologically viable model proposed in Ref. 



see also Ref. 



40|) 



f{R) = R-XoRo 1- 



i?2 



in which XQ,RQ,n are positive constants, satisfies the condition (11 01) . For this model, we 
have /'(O) = 1, i.e. coincidence with general relativity, again. However, because in this 
model /"(O) < 0, it is unstable |41|- 



The model described by 



f{R) = R + -ln{cosh{aR) + bsmh{aR)), 



where a, b are constants, admits the solution, too. Here, we have 

1 1 



/'(o) i + b- 



An example of the models incompatible with the condition (1T0|) is the so-called IR- 



modified gravity model of Refs. 



43 



44| described by 



f{R) = R- 



R 



where ^ ~ i/o, Hq being the Hubble constant. However, such models are ruled out by the 
Dolgov-Kawasaki instability 20||. It is interesting to note that in fact the above model and 
the Starobinsky model (the second example model discussed above) are special cases of a 
more general model given by 

/(i?) = i?-(l-nV (^^ 



which have been investigated in Ref. 45 1 



Another model which does not admit the solution is given by 46 1 



f{R) = R + a\n{^ ] + ^R"', 



a, /?, m being constants. 



V. GRAVITATIONAL WAVES IN NONMINIMAL J='{g) GRAVITY 



In this section, we consider a nonminimal J^{G) gravity model described by the following 
action [l^ 

S = j ^g{^R + J^iG) + (1 + K'H{g))Cn?j d^x, (14) 
in which an arbitrary function of the Gauss-Bonnet invariant 



g = R^- AR^.R'^" + R^^xaR""^" 



(15) 



is coupled to the matter field Lagrangian Cm with a coupling constant k. This represents a 
generalized version of the actions introduced in Refs. UJ, Il2 |. 

The equation of motion resulting from the above action is given by 



(16) 



where 



dg ' 



and 



-2g^"' R^'N p + 2R'"'VoN^ 

Now, we consider the pp-wave anstaz, Eq. ([1]). First, we note that the Gauss-Bonnet 
invariant, Eq. (ITSl) above, vanishes identically for the spacetime described by Eq. ([1]). For 
vacuum, the above field equation reduces to 

^ + ^) ^' ~ 2-^(0)^^^' = 0' (1^) 

and 

J'(O) = 0. (18) 

Thus, the vacuum solution is admitted only if the above condition holds. Then, Eq. f ITT]) 
reduces to the same equation governing the waveform in general relativity. Hence, the 
waveform is the same as the general relativistic counterpart. 

Back to the nonminimal coupling to a massless particle, by taking Cp as the matter 
Lagrangian and inserting the associated energy-momentum tensor into the field equation 
( !T6|) . we obtain the plane wave solution ([T]) whenever Eq. (fTSl) is satisfied. The relevant 
waveform is given by 

Kg{u, X, y) = il + kH{Q))Kas{u, x, y). (19) 

If, in addition, the function l-i{G) is chosen so that it satisfies 'H(O) = 0, then the waveform 
is not distinguishable from the general relativistic one. 



VI. GRAVITATIONAL WAVES IN F{G) GRAVITY 



The action for J^{G) gravity is given by 



S 



d'^Xy/^ 



(20) 



which is a particular case of the nonminimal model discussed in the previous section with 
K = 0. Thus, all such models satisfying Eq. (fT8|) admit both the vacuum pp-wave solution 
and the AS solution. Now, from Eq. f|T9|) . with vanishing n, it is obvious that the waveform 
Kg{u, X, y) is the same as the one in general relativity. 
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Examples of the J-'{Q) models satisfying the requirement Eq. (1T8|) include 



with n > which is shown in Ref. 



47l | that it could also pass solar system tests for n < 0.074. 
Also, the following cosmologically viable models 

-a^iVg., (21) 

(22) 



J^{g) = tan ^ (^^^ - afiVg*, 
J^{Q) = /Ua/^ In cosh ~ tt/Uv^*, 



(23) 



proposed in Ref. [48[ admit the plane wave solution for a = 0. Here, /i and Q^, are positive 
constants. 

The model presented in Ref. 



491 provides another example which admits the wave 



solution. It is described by 



Fo + F,{{g-Goy^+' + Gl^+'y 

where Fq, Fi, Gq are constants. 

It is interesting to note that both f{R) and J^iG) theories considered above can be 
obtained from a more general theory described by the following action 50 1 



S 



-g{F{R,g) + Cm)d^x, 



(24) 



which can be seen from the associated field equation which admits the above plane wave 
solution when similar conditions to those introduced above are satisfied. 



VII. DISCUSSION AND CONCLUSIONS 



In this work, we studied plane-fronted gravitational waves with parallel rays in the context 
of some extended theories of gravity. We considered f{R) and modified Gauss-Bonnet gravity 
with minimal and nonminimal couplings to matter and showed that they admit the vacuum 
pp-wave solution and also an Aichelburg-Sexl-type solution if certain conditions are satisfied. 
For f{R) gravity, the required condition is that f{R) vanishes for vanishing scalar curvature. 
This condition is compatible with the requirements for such theories to pass local gravity 
tests. It is also the same condition for a given model to admit the Schwarzschild solution. 
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Thus, for those models admitting the Schwarzschild solution, it should be possible to obtain 
the AS plane wave solution by boosting the black hole one as in general relativity. The 
explicit form of the wave profile depends on /'(O) and coincides with the general relativistic 
wave profile for some specific models including the well-known Starobinsky model. A similar 
condition holds for the modified Gauss-Bonnet gravity. In the latter case, the solution is the 
same as the general relativistic one. This was examined for several specific cosmologically 
viable models. For models with nonminimal coupling between the matter and geometry, 
more interesting options are available, including the possibility of (dis) appearance of (the) a 
contribution from the nonminimal coupling by choosing appropriate coupled function. 

The gravitational wave solution presented here might be used as an experimentally 
testable, at least in principle, criterion to distinguish between various extended gravity 
models. This can be achieved by looking at the behavior of two pointhke objects in the 
gravitational field of the massless source. For models with minimal coupling, this can be 
seen by measuring the relative acceleration of two nearby test particles separated by 
which is obtained from the geodesic deviation equation 

= -R^^.^x'^n-'x^ (25) 

and noting that for the spacetime under consideration, the components of the Riemann 
curvature tensor are proportional to the second derivatives of the wave profile K{u, x, y) 
with respect to the transverse coordinates. For models with nonminimal couplings where the 
particles do not move along geodesies as a result of extra forces coming from the coupling, the 
above equation should be modified by adding the relevant terms. However, since the extra 
force is proportional to the gradient of the scalar curvature or the Gauss-Bonnet invariant, 
both of which vanish in the above considered spacetime, the same argument still holds. 
The fact that the waveform obtained by application of the extended theories of gravity is 
proportional to the general relativistic waveform would then reflects in the observed relative 
accelerations. 

The results obtained here might also be used to study the scattering of particles at high 
energies in the framework of extended theories of gravity. Possible interesting extensions of 
the present work include a study of shock waves due to massless particles moving in curved 
backgrounds, and particles with arbitrary multipoles moving in curved spacetimes in the 
presence of matter fields and-/or a cosmological constant. 
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